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Effects of impurity scatterings on the strong coupling (SC) contribution, stabilizing the ABM 
(axial) pairing state, to the quartic term of the Ginzburg-Landau (GL) free energy of superfluid 
. . . '^He are theoretically studied to examine recent observations suggestive of an anomalously small SC 

' effect in superfluid ^He in aerogels. To study the SC corrections, two approaches are used. One is 

\ based on a perturbation in the short-range repulsive interaction, and the other is a phenomenological 

^ ■ approach used previously for the bulk liquid by Sauls and Serene [Phys.Rev.B 24, 183 (1981)]. It is 

found that the impurity scattering favors the BW pairing state and shrinks the region of the ABM 
pairing state in the T-P phase diagram. In the phenomenological approach, the resulting shrinkage 
'~{ , of the ABM region is especially substantial and, if assuming an anisotropy over a large scale in 

aerogel, leads to justifying the phase diagrams determined experimentally. 
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PACS numbers: 67.57.Bc, 67.57.Pq 



I. INTRODUCTION 



o . 

^ . The aerogel is a highly porous medium composed of Si02 strands with a porosity of about 98% in the case of 
5-H ' samples used in many experiments on superfluid "^He. A structural correlation length, ^a, of aerogel is identified 
with the typical distance between strands, i.e., ~ 30 — lOOnm. In the system of superfluid ^He in aerogel, the 
^ ' coherence length at zero temperature competes with in magnitude, and it will be necessary to consider two 
^ ; \ different regimes separately to understand roles of aerogel in this system. At shorter scales than ^a, the aerogel acts 
• as an anisotropic scatterer inducing an uniaxially anisotropic scattering amplitude of '^He quasiparticles, while it acts 
' as point scatterers for them at much larger scales than In any case, the aerogel plays roles of a pair-breaker for 
I . superfluid '^He with a p-wave pairing. 
■ O ' In liquid ^He in aerogel, it is now believed that (at least) two superfluid states occur. One is the A-like phase, 
^ , and the other is the B-like phase P, The latter has essentially the same properties as the bulk B phase and 
thus, is believed to be in the BW pairing state. In contrast, the A-like phase with an equal-spin pairing (ESP) has 
remarkably different features in the phase diagram from those of the bulk A phase which is in the ABM (axial) 
pairing state. Among them, remarkable ones are 1) a shift of polycritical point (PCP) to lower pressures, 2) the 
positive slope {dTAB/dP > 0), appearing even at high pressures, of the A-B transition curve Tab(^) in contrast to 
^ the negative slope in the bulk case [1, U, and 3) an anomalously high A1-A2 transition field Based on the 
(y-j feature 1), a different pairing state from the ABM one was proposed [7]. However, it has been clarified recently that 
Qs^ , 1) can be understood, under the assumption that the A-like phase is in the ABM pairing state, as a consequence of 
' an anisotropy existing over long distances in the aerogel sample [sj or of an act of aerogel as a quenched disorder [§] . 

I , Although the features at lower pressures on the phase diagram in aerogel are well understood within these approaches 

' [1, the features 2) and 3) pecuhar to the high pressure region remain to be explained. However, the point to be 
f — . ' searched for seems to be clear: It is well understood [l^ [ll| that the ABM pairing state in bulk is stabilized by a 
correction term to the quartic term in the Ginzburg-Landau (GL) free energy induced by the repulsive interaction 
in the normal state. This correction term, called as the strong-coupling (SC) correction, is enhanced with increasing 
the pressure P, because the effective repulsive interaction between quasiparticles will be stronger at higher P-values 
where the bulk A phase is present over a wider temperature range. On the other hand, in the liquid '^He in aerogel, 
the SC correction was estimated based on the observation 3) to be anomalously small, and such a small SC correction 
may be consistent with the positive slope in 2), because the positive slope may be realized if the temperaure width 
of the A-like phase at the high pressure end of the liquid phase is narrow enough [1, So far, the SC correction in 
the impure superfiuid '^He was incorporated in theoretical calculations using a simplified treatment in which impurity 
effects have been incorporated merely through a relaxation rate of quasiparticles J^. However, the observations 2) 
and 3) seem to imply that this relaxation time approximation is insufficient if the A-like phase is in the ABM pairing 
state. 

In this paper, the impurity effects on the SC correction are studied in detail and carefully on the basis of two analyt- 
ical methods, a purely diagrammatic approach [isj treating the four-point vertex (FPV) of the normal quasiparticles 
perturbatively in the strength of a bare repulsive interaction, and a phenomenological approach [141] in which the FPV 
is parametrized in terms of Landau parameters estimated from properties in the normal liquid ^He. We find that, in 
both of the two approaches, the SC correction is reduced by incorporating an impurity-induced new process in the 
SC correction, and the resulting shrinkage of the region of the ABM pairing state in the impure case is much more 
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remarkable in the phenomenological approach, suggesting that this approach [lj| is more suitable to a description of 
real liquid "^He. 

In sec. II, the details of our theoretical analysis are explained, and the resulting general formulation is applied in the 
ensuing two subsections to examine the phase diagram of weakly disordered '^He according to the two approaches. In 
sec. Ill, the obtained results are used to discuss the experimentally determined phase diagrams of the liquid '^He in aero- 
gel by incorporating possible anisotropics in the impurity scattering provided by the aerogel structures. Concluding 
remarks are given in the final section. 



II. MICROSCOPIC CALCULATION OF STRONG COUPLING CORRECTION 

The superfluid '^He is in a spin triplet p-wave pairing state with a gap function of the form 

Aa/3(P) = A(p) • {iaa2)al3 = A^.jPj {i<y 1^0-2) a/3, (1) 

where a = (tii, f72, fa) are Pauli matrices, p is the unit vector p/pp parallel to the momentum p, and pp is the Fermi 
momentum. Since the aerogel structure acts as a nonmagnetic impurity in superfluid '^He, we start from the BCS 
Hamiltonian with an attractive p-wave pairing interaction and an impurity scattering term, i.e., Hbcs = Hp + -ffimp- 
Here, 
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Oj-.^lq) = X! ~"-P+q/2."(^^A'C^2)Q/3ap+q/2,/3, (2) 

where ctp^a is the Fourier transform of the field operator ^'^(r) of quasiparticles, u{r) is an impurity scattering 
potential with a Gaussian ensemble defined by u{r) = and u(ri)M(r2) = (27rriV(0))~"'^(5(ri — r2), the overbar implies 
the random average, r is the life time of a quasiparticle, and iV(0) is the density of states per spin on the Fermi 
surface. The p-wave pairing channel has been assumed to be dominant in the interaction Hamiltonian 

Hqp= f I [V^'\v-v')n{v)n{v') + V^^\v-v')s^.{v)s^{v')], (3) 



between quasiparticles conserving the total spin, where n and are the density and spin-density operators of quasi- 
particles, respectively. 

Throughout this paper, we work in the mean field approximation where the roles [9] of the impurity scattering 
acting as a quenched disorder on the superfluid order parameter Apj (r) are neglected, and thus, spatial variations of 
may be neglected in considering the equilibrium properties. Further, in examining impurity effects on the SC 
correction, a possible anisotropy in the scattering events will be neglected since, as well as the anisotropy, the SC 
correction itself is a small contribution to the condensation energy. Inclusion of anisotropy of scattering events will 
be postponed to sec. HI in which results comparable with an experimental phase diagram will be shown. Then, we 
only have to examine the coefficients of GL free energy functional 

^GL - ^ (aA;^,A,,,, +/3i|A;,,A^,,|2 -(-/32(A;^,A^,0' +/33A;,^^:,,^MJ^-J 
microscopically. When the Born approximation for u{v) is used, and the s-wave component of the scattering amplitude 
is assumed to be dominant, the coefficient a in eq.(4) may be well approximated by the familiar result 

« = i^(0)[ln|;^+^(i + ^)-^(i)]. 



where Tco is the superfluid transition temperature of the bulk liquid in the mean field approximation. Note that, 
due to the p-dependence carried by the the pair-field vertex, the impurity-induced vertex sketched in Fig. 1(b) can be 
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FIG. 1: (a) Vertex (shaded part) between two quasiparticles and the pair-field A. (b) Example of the vertex with an impurity- 
induced correction which vanishes in the p-wave pairing when the scattering amplitude is independent of the momenta. 
The dashed line with a cross carries {2-kN{0)t)~^ . 



neglected, where a solid line is the quasiparticle Green's function and a dashed line with a cross symbol represents 
the impurity scattering. In the weak coupling approximation, other coefficients in eq.(4), obtained consistently with 
the above a, are given by 

/5o(r),,,(2).i , 1 



7C(3) ^2 AtttT'' 
7C(3)jV(0) 

MT) = 240..^. , (5) 

where (z) is diagamma function. 

The most stable pairing state is determined by comparing the condensation energy f'^°'^'^ — _q,2^(4^j^) between 
different pairing states with one another, where the index N indicates each pairing state. Since a is common to different 
pairing states, the stability depends only on /3n- We investigate the stability of the ABM pairing state relative to the 
BW one, that is, the relative difference between Pabm = /3245 and /3bw = P12 + /3345/3, where /3i2 = /3i + /32, and 

First, we consider the SC correction to Pj in clean limit based on the general framework given by Rainer and Serene 
(RS) [ig. We will not use the conventional spin fluctuation model [ll[ , because this model does not explain a property 
of the bulk B-phase: As can be seen in Table 5.1 of Ref . [TH . the pressure dependence of (3bw in the spin fluctuation 
model is opposite to the result [l4 [based on the RS's method and on the use of experimentally determined Landau 
parameters. The latter approach 14 1 seems to have explained the pressure dependence of the speciflc heat jump at 
the normal to B-phase transition 16 1. According to Ref. pJOj . the SC contribution to Pj, 6(3j, can be obtained up to 
the lowest order in T / Ep by examining diagrams raised in Fig[21 where an open square is the quasiparticle 4-point 



vetex (FPV) function in the T = limit, F^^''] r(pi,p2;p3,p4). Here, the renormalzed interaction Hamiltonian was 



assumed to take the form llj 



H'iT^ = \ I I I I *l(ri)*^(r2)rX.(ri,r2;r4,r3)vI/,(r4)vI/,(r3). (6) 

a,^,7,5'^''i "^""2 "'''3 "'""4 

Among the diagrams in FigO the former two diagrams can be neglected : Fig. 2(A) belongs to the weak coupling 
diagram and expresses an interaction-induced rcnormalization of a pair-fleld verte x ap pearing commonly to all super- 
fluid pairing states. This vertex correction has been examined in eq.(3.7) of Ref. [10] and estimated to be negligibly 
small compared with the SC corrections. Further, Fig. 2(B) vanishes after carrying out frequency summations. Then, 
we have only to calculate the remaining four diagrams. Their contributions to the GL free energy functional are 
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FIG. 2: Diagrams describing SC corrections in clean limit. Here, any frequency dependence in the FPV (the open square) is 
neglected. Diagrams (A) and (B) can be neglected for the reasons explained in the text. 



expressed as 

Fig.(C) = -ir^ E / / / G(^^(Pi)GL°i,(~Pi)G(°)(p,)GLt(-P2)G(?(P3)GLt(-P3)G(?(^^ 

ni,n2,n3 "'Pi "'P2 "'P3 

xri'^ila2;a3,a4(Pl' P2; P3, P4)r^''^^,/3, ./33,^, ("Pi , -P2 ; "PS, -p4)Aai/3i (pi ) Aa^fe (P2 ) A)33„3 (ps) A|34„ Jj54) , 

Fig.(D) = E / / / G(°)(pOGL°i,(-pi)[G(°)(p.)]2GLl(-P2)G(°)(p3)GLl(-P3)G(°)(p4) 

ni,n2,n3 Pi P^ P-' 

><ri'*ila2;a3,a4(Pl'P2;P3,P4)rg^,„,;^3,^,(-Pl,P4;-P3,P2)A„,^,(pi)A„,^(p2)Aj^^^^ 

Fig.(E) = -\t^ I I I [Gi°)(Pi)]^GL°],(-pi)[Gg)(p3)]^GLl(-p3)Gg)(p,)G(°)(p4) 

ni,n2,n3 Pi P^ "^P^ 

xri'^ila2;a3,a4(Pl' P2; P4, P3)r{j'^^,a3;a2 ,/34 ' P4; P2, Pi ) Aai7 (Pl ) A^^^ Jpi ) A^;,5 (^3) AJ^^ (^3) , 

Fig.(F) = -\t^ I I I [Gi°HPi)]^GL°i,(-Pr)[G(^)(P2)]^GLl(-P2)G(^)(p3)G(°)(p4) 

"i,n2,n3 -^Pi -^Ps -^Pa 

Xri'^?,a2;a3,a4(Pl' P2; P3, P4)ri''3,a4;/3i,ft(P3' P4; Pi : P2 ) Aai7 (pi ) A^^^^ {Pl) ^a^siPi) ^\p^{Pi) , 

where G^e}{v) is the quasiparticle Green's function in the normsal state in clean limit, P4 = Pi + P2 — P3, and the 
integral means J(fip/{2Tr)^. 

Now, let us turn to the case with effects of impurity scattering. In each process accompanying impurity scatterings, 
we will keep only the lowest order term in the impurity scattering: In fact, as a consequence of the process dipicted 
in Fig. 1(b), the non s-wave paired superfluid is destroyed when tTcq is of order unity so that we need to focus on 
the case of a weak impurity scattering in which tTcq > 1. Further, the assumption of a weak impurity scattering 
is reasonable in the sense that the aerogel is a porous medium with so dilute scattering centers. First, for single 
particle properties such as the selfenergy term of the quasiparticle Green's function, we use the Born approximation 
by neglecting multiple scattering processes. Then, this self energy correction can be incorporated in the expressions of 
the GL terms by replacing the Matsubara frequency e„ in Gi';^^(p) = G(°^(£„) by e„ = £„( 1 + l/(2r|£:„|) ). Hereafter, 
a quasiparticle Green's function will be expressed by the resulting one G^^ = G(°^(e„). If other impurity effects are 
not considered, the resulting expressions of SPj in this relaxation time approximation are equivalent to those used in 
previous works [l^ • 

Here, let us first discuss mean field phase diagrams following from this relaxation time approximation. When 
determining an AB transition curve Tab{P)^ the next 6th order term has to be taken into account in the GL free 
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energy in the weak coupling approximation where different coefficients PjS have the same temperature dependence 
proportional to as one another. However, the SC correction to (3j brings additional dependences on T and P, 
and, at least close to Tc{P), a r4 R(P )-curve may be determined by eq.(4) with the SC correction but with no 6th 
order term. As in previous works TT], 1^ determining the bulk Tab line, it is natural to expect both the 6th order 
term and the SC correction to have to be incorporated in calculations. However, we have found that, even for the 
bulk ■^He with Tab{P) far apart from Tco{P), reasonable Tab{P) curves can be obtained from eq.(4) with no 6th 
order term in the two approaches to be explained later. For this reason, we determine hereafter the phase diagram 
according to the GL free energy, eq.(4), truncated at the 4th order in A^,i. 

In Fig. 3, dependences of the width of the A-like region on {tTcq)"^ resulting from this relaxation time approximation 
are shown. There, the Tab{P) curve is found to rather decrease with increasing the disorder strength {tTco)~^ ■ It 
implies that the SC correction in this relaxation time approximation is rather enhanced with increasing disorder: 
Based on eq.(4), the Tab line is determined by the relation /3abm = /SbWj with the SC corrections to f3j included, 
which is independent of Tc{P) determined from a in eq.(4). Hence, the impurity-induced reduction of Tab seen in 
the r-dependence of Fig. 3 implies an extension of the ABM state region if the familiar impurity-induced reduction of 
Tc{P) is absent. Once considering correspondences with experimental facts suggesting a remarkable impurity- induced 
reduction of the SC parameter [Q| (see sec. I), we feel that an impurity- induced reduction of the SC parameter will be 
present and needed to describe superfluid '^He in aerogel properly. 
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FIG. 3: Mean field phase diagrams in the case with no anisotropy and in the relaxation time approximation obtained in terms 
of (27rr)"^ = 0.08 (left), 0.13, and 0.19 (right) (mK). The shaded regions denote those of the mean field A-like phase appearing 
without efi'ects of anisotropies and the quenched disorder. The AB transition curve afi'ected by the impurity scatterings (solid 
curve) is shifted to /ouier temperatures as the impurity concentration ~ {Tcot)~^ is increased. For comparison, the bulk AB 
transition line determined experimentally is indicated by a dashed curve in each case. 



For the reason mentioned above, we investigate hereafter impurity-induced corrections to the two-particle processes, 
which include impurity- induced vertex corrections and will be denoted as 6Pj later. Among such vertex correction 
terms, the lowest order terms in (tTco)^^ contributing to Sf3j are obtained by replacing one of two FPVs in Fig.2C - 
2F with a bare impurity vertex or adding an impurity line to Fig. 2 A and 2B. The resulting nonvanishing diagrams 
are those in Figd) All diagrams in Fig|3]are obtained by adding a single impurity line to Fig.2B. On the other hand, 
the corresponding ones occurring from Fig.2A, which are sketched in Fig. 5, vanish or can be neglected: Figure 5(a) 
can be absorbed into Fig.2A itself which, as already mentioned, is negligible. The next Fig. 5(b) is absorbed into an 
additional weak-coupling quartic term including a single impurity line which was denoted as — £imp/3o in Ref.0]. As 
noted there [^, however, this ffimp-term does not contribute at all to the relative stability between different pairing 
states and hence, can be neglected in determining Tab{P) according to eg .(4). Further, Fig. 5(c) is found to vanish at 
least up to the lowest order in T/Ep- This is verified by, following Ref.[l3|, replacing the Green's function and their 
products appearing outside the FPVs in the way 
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Gap)G_e(-p) 
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and performing necessary momentum integrals and frequency summations. Further, according to the "rule" of the 
T/i?i?-expansion formulated in Ref . [lol| . the corresponding diagrams arising from Fig.2C-F and with a single impurity 
line are of higher order in T / Ep compared with those of Fig. 4 because an impurity line plays roles of a FPV with no 
finite frequencies carried. In this way, for our purposes of examining the SC corrections to (jj up to the lowest order 
in T/Ep: we only have to focus on the diagrams in Fig. 4 as the impurity effects on the SC correction which were not 
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FIG. 4: Diagrams describing the impurity-induced SC corrections which do not appear in the relaxation time approximation. 
Note that the diagrams (a), (b), and (d) are identical with (i), (f) and (h), respectively. 

considered in the relaxation time approximation. Their contributions to the quartic term of the GL free energy are 
given by 



Fig(a) 



Fig(b) = 



I^^E^E / / / G2^(Pl)G-e,(-pl)G,,(p3)G^Jp4)G_,,(-p4)G,,(p2) 
ni n2 "^Pi "^P2 "^P3 

X |W(P1 - P3)\'^T'^ala2:a,,ai (Pl, ^1, P2, £2; PS, £l, P4, 62) A^^ffipi) Aj^^^ (pi) A^^jipi) A\^^ (^4), 

\tY,TY.I I [ Ge,(pi)G_,,(-pi)G,,(p3)G_e,(-p3)G2^(p4)G_e,(-p4)Ge,(p2) 
m 712 -^Pi -^P^ "'P3 

xKPi-P3)prW (pi, £1, P2, £2; P3, £1, P4, £2)A„,^^(pi)A|^„.^(p3)Aa,^(p4)Ajy„Jp4), 



Fig(c) = ^tj^tj: 



G^^ (p3)G_e, (-P3)G,, (pi)G2^ (P4)G_,, (-P4)G,, (P2) 



Pi •'P2 •'P3 



X \u 



Fig(d) 



Fig(e) 



Fig(g) = 



(Pi-P3)prw 

,a2;a3,Q:4 

(pi, £1, P2, £2; P3, £1, P4, £2)Aa,/3(j53)A|,„.^(p3)Aa27(p4)At„^(p4), 

ir^T^ / / / G2^(pi)G_,,(-pi)G,,(p3)G,,(p4)G_,,(-p4)G,,(p2)G_,,(-p2) 

ni n2 "'Pi "'P^ ■'Ps 

xKpi-p3)|^r(f) 

,a2;Q:3,a4 

(pi,£l,P2,£2;P3,£l,P4, £2)A„,^(pi)A^„3(pi)A„,^(p2)At^Jp4), 
\tY.TY,I I I G,,(pi)G_e,(-pi)Ge,(p3)G_e,(-p3)G,,(p2)G_,,(-p2)Ge,(p4)G_e,(-p4) 

X |u(pi - P3)Pr(f^)_„^.„3_„^ (Pl, £1, P2, £2; P3, £1, P4, £2)A„i0(pi) A^^3 (p3)A„,^(p2) AJ^^^ (^4), 

\tY,TY. Ill G?^(pi)G_,,(-pi)G,,(p3)G2^(p2)G_e,(-p2)Ge,(p4) 

m na "'Pi "'P2 "'Ps 

(Pi-P3)prw 

,0:2 ;a3,a4 

(Pi, ei, P2, £2; Pa, ei, p4, e2)Aai;3(pi)A^^3 (j5i)Ac<27(p2)Ajy„Jp2). 



X \u 



(a) (b) (c) 

FIG. 5: Impurity-induced diagrams occurring from Fig.2A which can be neglected in the present analysis. 



Hereafter, let us construct formulae of two terms, SPj and 6(3j, giving impurity-dependent SC corrections to (3j 
which are useful numerically once an appropriate FPV is substituted. Further, as far as the lowest order contributions 
in T/Ep are concerned, one can neglect any frequency dependence of FPVs in Fig. 4, and the integral over 

= P^/{2m) — 11 can be simplified using eq.(7). Next, the renormalized FPV, i.e., the scattering amplitude due to 
mutual interactions between the quasiparticles, will be decomposed fl?} into a spin-symmetric part T^*^ and a 
spin-antisymmetric part 

2^(0)r53;^, = T'^'^ 5^,5 PS + T^^Hc^U ■ i'f)0S- (8) 
Alternatively, it can also be expressed by the spin singlet amplitude and the triplet amplitude Tt as 

where 

Further, by assuming the bare impurity scattering amplitude |u(p)p to be dominated by the s-wave component and 
separating the momentum integrals from the frequency summations, the diagrams in Fig. 2 and Fig. 4 are rewritten as 

^. N{0)~^ T fdn, rd^2 [d^. 

x([T(^)(pi,P2;P3,P4)]'{[A*(pi) • A(p3)][A*(p2) • K{p\)] 

-[A(p3) • A(p4)][A*(j5i) • A*(p2)] + [A*(p2) • A(p3)][A*(pi) • K{p^)]] 

+ [r(^)(pi,P2;P3,P4)]'{ - 5[A*(pi) • A(p3)][A*(p2) • A(p4)] 

+5[A(p3) • A(p4)][A*(pi) • A*(p2)] + 3[A*(p2) • HpM^*{pi) ■ Mpi)]}), 

X ([TW(pi,P2;P3,P4)T(^nP3, ~P2;Pi,-Pi) 

+T^''\p\,P2;P3,P4)T^''\p3, -P2;P1,~P4)] |A(pi)HA*(p2) • A(p4)] 

+ [t'^''\pi,P2;p3,P4:)t'^''\p3, -P2;pi, -m) + t'^'^\pi,P2;p3,P4:)t''''\p3, -P2;pi, -Pi)] 
x{[A(pi) • A*(p2)][A*(pi) • A(p4)] - [A(pi) • A(p4)][A*(pi) • A*(p2)]}), 

x[[T'^^\p\,pr,P3,p\)]\\A{pi)\'\A{p3)\' - |A(pi) . A(p3)|' + \A{n) ■ A*{p,)\^] 
+ [r(^'(Pi,P2;P3,P4)]'[3|A(pi)|2|A(p3)P + \A{Pi) ■ A{P3)\' - |A(pi) • A*(p3)P]), 
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Fig(F) = - ^5.- j—j^j ^K\P. -PsI - 1) 

X ([T(^)(pi,p2:P3,P4)]'[|A(pi)nA(p2)|2 - |A(pi) • A(p2)P + \^{pi) ■ ^*{P2)?] 
+ [t(^Hpi,P2;P3,P4)]'[3|A(pi)|2|A(p2)P + |A(pi) . A(p2)|' - |A(pi) • A*(p2)|']), 

A ^(o)V v^ 1 \^ 1 rdfti fdn^ fdn^ 
= I^(S (n + 1 + ^)0 177 17/ 

x(T^^\p\,pr,P3,P4)\A{p,)\'\A{p^)\'' 

+T(-)(pi,P2;P3,P4){|A(pi) • A*(p4)|' - |A(pi) • A(p4)P}), 
X (t(«)(pi,p2;P3,P4) [A(pi) • A*(P3)] |A(P4)|' 

+T(-)(pi,P2;P3,P4){ [A(pi) • A*(p4)] [A(p4) • A*{ps)] - [A{p,) ■ A{pi)] [A*{p3) ■ A*{pi)] }) , 
x(T(^)(pi,P2;P3,P^4)|A(p3)nA(p4)|' 

+r(^)(pi,P2;P3,P4){|A(p3) • A*(p4)P - |A(p3) • A(p4)P}), 

F^g(d) =- 8;^(g (n+l + ^)0 7^/ 17/ 17/ 17^(b^+---|-^) 

x(tW(pi,P2;P3,P4)[A(p2) • A*(p4)]|A(pi)|2 

+T(^)(pi,P2;P3,P4){[A(pi) • A*(p4)] [A(p2) • A*(pi)] - [A(pi) • A(p2)] [A*(pi) • A*(p4)]}), 
^''^'^ = 4^(g (.+ i + ^)0 7^/17/17/17^(1-+---!-^) 

X (r(^)(pi,P2;P3,P4)[A(f)i) • A*(P3)] [A(P2) • A*(P4)] 

+T(-)(pi,P2;P3,P4){ [A(pi) • A*(P4)] [A(P2) • A*(P3)] - [A(pi) • A(p2)] [A*(p3) • A*(P4)] }) , 

N{<df (^ 1 ^2 1 fdn, fdn2 fdns 
F^s(g) = T^(i:(77IT3;j^) 7^^/17/17/17^(1-+---'-^) 

x(t(«)(P1,P2;P3,P4)|A(pi)|'|A(P2)|' 

+tW(Pi,P2;P3,P4){|A(pi) • A*(p2)|' - |A(pi) • A(p2)|'}), (11) 

where the factors such as Sa denote results of summations over Matsubara frequencies; 

Se = Sf = (TrfcsT)"^ V V V ^^^^^^^^[sgne3][sgn(ei +£2 - £3)], 
— — — £1 £2 

ei £2 £3 

^ \e\\ £2 £3 
£1 £2 £3 ' ' ' ' ' ' 

~Sc = {-ksTrT.i:i:ih\h\h \i J eV 
'-^ '-^ '-^ £1 £2 £3 £1 + £2 — £3 

£1 £2 £3 

The remainder of our calculation is to perform angle averages over the Fermi surface. To conveniently parame- 
terize T'^**)(p\,P2;P3,P4) and r'^")(p\,p2;P3,P4), which are functions of only two independent variables, we introduce 
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Abrikosov-Khalatnikov angles 6 and 4>, which are related to the angles between incident two particles or between an 
incident and an outgoing particles in the manner 

Pi ■ P2 = P3 ■ Pi = cos 6* EE Xi , 
2 . 2 ^ I 

Pi-Pz^ COS - + sni - cos = X2, 

2 ^ • 2 ^ / 

Pi ■ Pi = COS - — sm - cos = X3. 

To treat the angle-integrals in eq.(ll) including a 5 function, we choose pi + p2 as the polar axis for p^, and the 
azimuthal angle of will be defined by measuring it from the plane containing pi and p2- Then, dils = d cos 03^0, 
and we can perform the integral with respect to 9^. Further, the unit vector p will be chosen along pi -\- P2, and an 
additional azimuthal angle if) is chosen as an angle measured, for the new polar axis p, from the plane containing pi 
and p2 ■ Then, we have [lO| 



For instance, the product of scattering amplitudes T{pi,p2;p3,P4)T{p3, —p2;pi, —pi) appearing in the expression for 
Fig.(D) in eq. (fTT|l can be expressed as T{d, 4))T{6' , cj)'), where 9' and (j)' are related to 9 and (j) through the formula 

cos9' = cos(/) - cos^(0/2)(cos(/) + 1), 

3cos2(6i/2) - 1 - [cos2(6'/2) - l]cos(/) 



cos2(6'/2) + 1 + [cos2(e'/2) - 1] COS<j) 
For later convenience, the variables x'l, X2, and x'^ will also be defined as 

x[ = cos 9' = — xs. 



X2 = cos — + sm — cos = X2, 
I 2 ^' -2^' 1 1 

Xn = COS sm — cosffl — —XT. 

^ 2 2 ^ 



What we need to perform finally is to calculate the integral 

dflp f dip 



[A*(p,)-Afe)][A*fe)-A(pO]- (12) 

47r J Zn 

Since the angles 9 and <j) fix the relative position of the four vectors pi{i — 1 — 4), and the absolute positions of 
those vectors are determined by p and tp, the polar coodinates {9i, (pi) oipi are (^1, ^1) — {9/2, 0), (^2, (P2) = (^/2, tt), 
(^3j 03) = (p): and (6*4, (pi) = {9/2, (p + n) in the body-fixed frame. Each component of A(pi) can be expanded in 
terms of Z = 1 spherical harmonics as 

I 

^f^{Pi) = ^ B^,jnYlm{pi) 
m—~l 

with 1 = 1, where, by definition of A^j, the relations between A^^^j and Bf^,j{l = 1) 

_ 1 [8^ 

I/Stt, ■ a \ 



B,,o = (13) 



are satisfied. Then, each integral appearing in the expressions of eq.(ll) reduces to 

^^l"-^(P^)^^"^M)YlmAPk)yimAPl)- (14) 
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To perform these angle integrals, it is convenient to use the formula 

I 



where R is the rotation which maps the coordinate system (pxyPy^Pz) into the body coordinate system parameterized 
by the angles 6i and i^i , and D^^, (R) is the corresponding rotation matrix. Using the standard properties of rotation 
matrices together with this transformation, we can obtain 



2L+1 

L=0 



X {lmil'm2\Lmi + 'm2){l'm^lmi\L — mi — m2)^^l\0,(j)), (15) 

where 



Note that the Clebsch-Gordan coefficients in eq. (|15p require rrii + m2 + + TO4 — and m'l + •m'2 + m'^ + m\ — 0. 
Using these expressions, one can obtain the contributions to the quartic term in the GL free energy functional. The 
same procedure has been used in clean limit in Ref . [To| . 

Results of the SC correction to fij in the relaxation time approximation, given by (C)-(F) in Fig. 2, are given by 

= Spf + 6(3f+'' + Spf , (16) 

where 

for K = C, E + F, where Se+f ~ Se — Sf, and 



+VP{0, 0) [T(^) {0, 0)T('^) {9\ 4>') + r(") [9, dp)T'^'^ [0', 4>')]). (17) 



= -^7^/3o(r)(ft3/4)(wf (0, 0)[r(^) (0, 0)r(^)(0', 0') + t(")(0, 0)r(")(0', 0')] 



The weighting functions Wj^ and V^^ are given in Table. I, and {Z) denotes the angular average 
/„^d(cos(0/2))/f Zd0/(2^). 



TABLE I: Weighting functions for Spf. 





K = C 


K = D 


K = E + F 




-4x? + + 





-fx? -3x1 + 1 




20a;? - 5x1 - 5a;i 


3xiX3 — X2 


f xf + 3x1 - 1 




-23;? + ixl + ixl 


2x1X3 — 4X2 


ix? + xl + 1 




2x1 - 23x1 + 17xi 


— 3X1X3 + X2 


-f-xl-9xl + f 




%xl - 2xi - 2xi 


— 3X1X3 + X2 


3xi + 6x1 — 3 




-8x? + 2xi + 2xi 


— 3X1X3 + X2 


-3x? + 6x1-1 




-2x1 + 3xi + 3xi 


— 3X1X3 + X2 


ix?+xl + f 




2x? + 17x1 - 23xi 


2x1X3 — 4X2 


fx? + 11x1 -f 




2x1 - 3x1 - 3xi 





1 ^2 2 9 

2X1 - X2 - 2 




-lOx? + 15x1 + 15xi 


X1X3 + 3X2 


fx? + xl + f 



On the other hand, the corresponding results of the SC corrections 5[3i including the impurity-induced vertex 
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correction, which arise from (a)-(i) of Fig.4, are 



7C(3)' ' ''iTTTTEpy ^2 At^tT 

ICi^y^^ ' A-ktT Ef 

1 T / .,,^,l 1 



where the weighting functions wf^^ and V^^"^^ are shown in Table. II. 



TABLE II: Weighting functions for 6(3^ 



(k) 





k = 


a 


k = b = d 


k = c 


= g 


k = e 






















-3a;i 


+ 1 


— 3x'lX';s + X2 


■J 2 

-3xi 


+ 1 


— 4x? + xi + xi 




-2x1 


+ 4 


4X2 — 2X1X3 


-2xf 


+ 4 


-x\ + 4x2 - 4x1 




3^1- 


- 1 


3X1X3 — X2 


3x? - 


- 1 


-X? - xi + 4x1 


wt^ 


3^1- 


- 1 


3xiX3 — X2 


3x? - 


- 1 


4xf - xl - xl 




3a;i- 


- 1 


3xiX3 — X2 


3x? - 


- 1 


4xf — x'l — x'l 


wt^ 


3x1- 


- 1 


3X1X3 — X2 


3xf - 


- 1 


-x'i -xl+ 4x1 




-2x1 


+ 4 


— 2xiX3 + 4X2 


-2x1 


+ 4 


2 1 ^ 2 2 

-Xi + 4X2 - X3 


wt^ 




















-4 


-3 


— X1X3 — 3X2 


-X? 


-3 


2xf — 3x1 — 3x1 



Summarizing the above results, we reach the formula of 5(3j 

47r2 ^ 1 T 



which represents the new impurity- induced SC correction and should be added to 50j. Here, the values of W*°*^' 
and y.tot'^-i are shown in Table.III. 

TABLE III: Weighting functions appeared in eq.(18). 



- ixl - 4x1X3 

^^(total) _2a;2_2a;l+ 4X1X3 

y^Ctotai) 1^2 _ 9^2 ^ 4^^^^ 

^Jtotal) _9^2^ 1^.2 ^43.^^3 
_9^2^ 1^2^4^^^3 
1^2_ 9^2^4^^^3 

^^(total) _2a;2_ 2x1 +4X1X3 

^^(total) Q 

(total) 3 2 I 13 2 o 
V5 f Xi + X3 - 8X1X3 
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Our remaining task is to determine the FPV concretely. To do this we use two approaches; one is a perturbative 
approach starting from a bare repulsive interaction [l3| and the other is a phenomenological approach in which the 
FPV is determined in terms of Landau parameters estimated from experimental data of transport coefficients and 
others [13]. 



A. Perturbative approach 



In this approach [l3| , the attractive interaction between the quasiparticles is assumed to be induced by a short-range 
repulsive interaction between two bare particles, and hence, we start from the following two particle interaction term 



(19) 



where g is a positive coupling constant. We carry out the perturbative expansion in g and calculate the FPV function 



up to the second order in g. 



ri/3;^5(Pi. P2; P3, Pi + P2 - Ps) = [^g + 2,g^n(p2 - ps) + g^C(pi + pa) - .g^n(pi - Pa)] Sa^Sfjs 

- [^.9 + .9'C(Pi + P2) + 5'n(pi - P3)] (a)„^ • {a)ps: (20) 



where [1; 



n(q) = -iV(0) 



- ,„ 2pF + |q| 



In- 



4pF|q| 2pF-|q| 



is a polarization function, 



C7(q) = N[0) 



2 vi - qV4p^ + 1 



is a Cooper loop regularized at a large momentum. Comparing with Eq.® we can obtain T'^^'^ and T^**^ with the use 
of the Abrikosov-Khalatnikov angles 9 and in the manner 

T^''\e,(^) = A + A2{27r(2;3) + c(rri)-^(x2)}, 

r(-)(0,0) = -A-A2{c(2;i)+7r(x2)}, (21) 
where A — gN{0), and the functions tt{x) and c{x), corresponding to n(q) and C(q), respectively, are given by 



, , x/2(l + x) , V2 + VT 
Tr{x) — 1 ^ : — ,^ In — ;= 

c{x) = 2| 



— X 

VI - a; V2 - \/\-x 

2V2 V2 + VT^ 



(22) 



As one can see from the formula (fT7|) . the combination of T'^'^^T'^^^ {i^ j — s,a) always appears in calculating the 
SC corrections to Pj. Hence, we can include third-order corrections in A to the quartic term of the GL free energy 
functional. The products of two vertex parts can be expressed as 

[tW(0,0)]' = X' + 2X^27:{x3)+c{x,)-7r{x2)}, 

[r(-)(0,0)]' = A2 + 2A3{c(xi)-t-^(x2)}, 

T(^)(0, 0)r('^) {e\ 0') = -A2 - X^{27r{x3) + c{xi) - tt{x2) + c{x[) + 7r(4)}, 

T(")(0,(/.)T('')(0',0') = -X''-X^c{xi)+Trix2) + 27:{x's)+c{x[)~n{x'2)}. (23) 

Then, SPj, i.e., the SC contribution to (3j in the relaxation time approximation, is given by 



T 



6/3j = -Po{T)—X''[{mfSc + mf+''SE+F+mfSD) 
+ X{nfSc + nf+^'SE+F + rifSo) 



(24) 
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TABLE IV: Numerical values of mf and nf in eg. ((24)) . 



j 


my 




my 


ny 




nf 


1 


1.0943 


0.2344 


0.3128 


-4.0586 


-1.4992 


-0.1950 


2 


0.0002 


5.1601 


-3.1273 


-0.7179 


12.3408 


-12.0611 


3 





0.4689 








1.1323 





4 


-0.0002 


1.0948 


3.1273 


0.7180 


-4.8988 


12.0611 


5 


2.1890 


2.0326 


-2.5016 


7.8083 


3.0325 


-12.4511 



where the coefficients in eg. ([24)) are given in Table IV. 

As the limit of vanishing r"^ of eq.(22), the corresponding expression in clean limit becomes 

= -(3o {T) ^ X' (M, + AiV, ) , (25) 

where Mi = 76.1, iVi = -271.3, Ah = 7.2, N2 ^ -119.6, M3 = 6.4, N3 = 15.5, A'U = 48.4, N4, = 221.5, A/5 = 110.3, 
and 7V5 = 264.1. It is found that the use of these numerical values of coefficients leads to reproducing the numerical 
results in Ref . [l3| . For the ABM pairing state, the coefficients are given by Mabm = M245 = 165.9 and A'abm — 366.0, 
while those for the BW pairing state are A/bw = + A/345/3 = 138.3 and A^bw — —223.9. The third order terms 
are important for the stability of the ABM pairing state because the sign of Abw is negative contrary to that of 
Aabm- 

On the other hand, the result on eq. (|18p in this approach is given by 

Sp^ = -^^f] m— ^f^(i)(i + _J_)Y(p +XQ ) (26) 

7C(3)'^°^ UtttTEf\ ^2 4tttT'J ^ ^ ^ ' 

where Pi = -0.81, Qi = 1.4, P2 = 0, Q2 = 0.002, P3 = 0, Q3 = -0.66, P4 ^Q,Qa = -1-8, P5 = -1.6 and Q5 = 0.45. 
That is, Pabm = —1.6 and Qabm — —1.32 while Pbw — —1-3 and Qbw — 0.81. The third order term in both A and 
1/r has, in the BW pairing state, the opposite sign to that of the ABM pairing one. That is, the impurity effects in 
the BW state enhance the SC correction, while they rather reduce that in the ABM pairing state. 

A typical result of T-P phase diagram in this perturbative approach is shown in Figl6] where A was assumed, 
for simplicity, to be independent of pressure. In obtaining a phase diagram based on our calculations, we use 
hereafter the experimental P-dependence (isj of the bulk transition temperature Tcq{P) and the P-dependence of 
the Fermi energy given in Ref.[i^. The dashed lines are the transition curves in clean limit (i.e., of the bulk liquid). 
The resulting bulk TAB(P)-curve is slightly wavy because the neglect of P-dependence of A and the perturbative 
treatment in A may not be fully justified. The lower solid curve with negative dP/dT values is the TAs-line obtained 
in the relaxation time approximation, i.e., by neglecting 5f3j. On the other hand, the upper (inner) solid curve close 
to the corresponding (dashed) curve of the bulk liquid is the Tab curve resulting from the full expressions including 
5(3j. Here, we have chosen the value l/(27rr) = 0.13 (mK). From the figure, one can see that the impurity-induced 
SC correction, 5(3j, weakenes the stability of the ABM pairing state. 



B. Phenomenological approach 

In this approach [l^ . the quasiparticle FPV function r(^''(pi,p2;/%,P4) depends on two momentum transfers qi = 
\pi — P3I = 'i.kpyj (1 — X2)/2 and q2 = \pi — Pi] — ^kp^ (1 — a;3)/2. The singlet and triplet scattering amplitudes are 
approximated by Ts(x2-,xj) = Ws{x2) + Ws{x^) and Tt{x2TX^) = Wt(x2) — Wt(x^), which automatically satisfy the 
exchange antisymmetry conditions, Ts{x2,x^) = ^5(^3, 0:2) and Tt{x2.,x^) = — Tf (2:3, X2). We expand Wg and Wt in 
terms of the partial waves 

00 

Ws{x) = Y.^iPi^^)^ 

1=0 

oc 

Wt{x) = Y^WlPiix), (27) 

1=0 
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Temperature(mK) 



FIG. 6; Mean field phase diagram of the model with no anisotropy obtained in the perturbative approach and in terms of 
l/27rr = 0.13 (mK) and A — 0.821. The dashed lines are the transition curves of bulk ^He. By including the impurity-induced 
contribution, ea. (|18|l . the stability of the ABM pairing state is weakened and, as indicated by the arrow, its region shrinks until 
reaching the shaded region. 



where Pi{x) is a Legendre polynomial, and determine the coefficients Wi'^ from experimental data. Note that Wq 
can be ignored in considering the physical quantities since / = component disappears from Tt{x2, xs). Following 
Ref. we use the quasiparticle lifetime, the thermal conductivity, the viscosity, and the spin diffusion to determine 
three Landau parameters Fq, Fq, and Ff through experimental estimations for them and to obtain the coefficients 

(1) Quasiparticle lifetime 

riO)= 4^^^^ ' 



7r3TO*fc|T2 {W{e,<j)))' 



(2) Thermal conductivity 



(3) Viscosity 



(4) Spin diffusion 



K = {7ry2)nkB{T/TF)vlT{0)SE{X.), 

A„ = {w{e,<j)){i + 2cose))/{w{e,(b)), 



■q = -nvFPFT{0)SoiXn), 



A,; = {W{9, - 3sm\9/2)sm'' ^))/{W{e,cj))), 



D = ^vUl + FS)r{0)So{\D), 



1-Xd = {W^i{0,^)sm^ie/2){l-cos^))/{Wie,ct>)), 
The scattering rates, W{9,(j)) and Wn{9,(j)), and the functions Se{o)W given by [l4 



Se{o) 



)2 + iT,(0,0)2 + irt(0,0)r,(0,0), 
)2 + iT,(0,0)2 + irt(0,0)r,(0,0), 



E 



2n+l 



n— cvcn(odd) 



n(n + l)[n{n + 1) - 2A] ' 



(5) Landau parameters Fq, Fq, and F^ 



l + Ff/i2l + l) 
Ft 

l + Ff/(2? + l) 
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where Wm(1) = Ez=oW^*''- 

We expand the scattering amphtude up to Z = 3, 

Ts{x2,x^) = 2W^ + Wt{Pi{x2) + Pi(a;3)) + W^{P2{x2) + ^2(2:3)) + W^g (^3(x2) + ^3(2^3)) 

rt(x2,X3) = Wl{Pi{x2) ~ Pl(x3)) + W^2*(^2(a;2) - P2{xz)) + WI{P^{X2) - P3{X3)). (28) 

Then, we have seven physical parameters, t{0)T^, A^, A^, Xd, Fq, Fq, and , and adjust the seven fitting parameters 
Wl''\l < 3) to minimize the sum of squared deviations of the calculated physical quantities from their corresponding 
experimental values. Using eq.(lO), we can express T^''-' and T^^'' in terms of Tt and Ts as 

T'^-'>{9,(t>) - i(rt(x2,X3)-T,(x2,X3)), T<^'He,<P) = \{m{x2,X3)+T,{x2,X3)), 

T^-\e',c^') = i(r*(4,4)-T,(x^,4)), T(^\e'A') = i(3r,(4,xi) + r,(4,4)), (29) 



and calculate the SC contribution 6(ij concretely. Including a frequency cutofT Ec in the gap function A in the form 
[1^ A(£„;T) = A(T)/[1 + {sn/sc)'^], we obtain the following expression of the SC correction to j3j in the relaxation 
time approximation 



5/3, 



T r 



E+F 



ifs 



D 



where the coefficients are given in Tables V 



(30) 



TABLE V: Numerical values of the coefficients / 



K 



P=12 (bar) 


P=16 (bar) 


P=20 (bar) 


j 


h 


■i 


'j 


h 


'.7 




h 


■i 


h 


1 


0.8543 


0.6598 


1.0695 


1.0746 


0.7405 


1.2171 


1.1026 


0.7588 


1.2137 


2 


-1.4707 


15.3879 


-2.4577 


-1.6729 


17.5695 


-2.8381 


-1.7178 


18.1001 


-3.0316 


3 


1.4660 


3.1556 


3.0504 


1.7383 


3.3996 


3.5672 


1.7339 


3.4850 


3.5957 


4 


3.9382 


1.7999 


6.8452 


4.5471 


1.9915 


8.0498 


4.6242 


1.9306 


8.2253 


5 


7.7364 


2.6753 


-1.6559 


8.2614 


2.8890 


-2.0485 


8.5912 


3.0719 


-2.2023 



P=24 (bar) 


P=28 (bar) 


P=34.4 (bar) 


j 


h 




h 


h 




h 


h 




h 


1 


1.2522 


0.9086 


1.0230 


1.3338 


0.9154 


0.9880 


1.1269 


0.9168 


0.9559 


2 


-2.1025 


19.9582 


-4.9303 


-2.1983 


21.3248 


-5.9904 


-2.1791 


20.0521 


-5.4358 


3 


1.1921 


3.3555 


2.8961 


0.9999 


3.3166 


2.5887 


0.8187 


3.3170 


2.3066 


4 


4.5375 


0.5847 


8.2523 


4.4206 


-0.2511 


8.6678 


4.2128 


-0.1348 


7.4739 


5 


10.2905 


4.5880 


-3.3103 


11.1117 


5.3985 


-4.1031 


11.1826 


5.2853 


-3.2555 



The SC correction to f3j in clean limit following from eq. (|30p is given by 

<5/3, = -/3o(T);^A/3,. (31) 

Numerical values of A/3j, in which the factor Tco/Ep is also included to collect the pressure dependence, are shown 
in Table fVll where we set the cutoff Sc = 0.068£'f. These numerical values are almost same result as that in Ref . [l3| . 
As one can see from the relative difference |A/3abm — A/3bw|, PCP is obtained to be 24 bar which is close to the 
experimental PCP, 22 bar. 

In the present case, the new impurity- induced SC correction to Pj, ea. (|18p . is given by 
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TABLE VI: Numerical results of A/3j in eq.(31). 



P[bar] 


A/3i A/32 


A/33 


A/34 


A/35 


A/3abm 


A/3bw 


12 


0.043 0.066 


0.105 


0.208 


0.229 


0.503 


0.289 


16 


0.057 0.087 0.135 


0.267 0.267 


0.621 


0.367 


20 


0.063 0.099 


0.147 0.293 


0.298 


0.690 


0.408 


24 


0.071 0.088 


0.125 


0.291 


0.376 


0.755 


0.423 


28 


0.077 0.092 


0.119 


0.295 


0.421 


0.808 


0.447 


34.4 


0.073 0.096 


0.115 


0.288 


0.461 


0.845 


0.457 



where A/?^- are shown in Table IVlTl In comparison with Table fVTl all signs of A(]j are opposite to those of A/3j , implying 
that the impurity scattering weakens the SC corrections for both the ABM and BW state. As for the relative stability 
of the pairing states, the impurity scattering favors the BW pairing state since the relation |A/?abm| > |A/3bw| is 
always satisfied. 



TABLE VII: Numerical results on A/3j in eq.(32). 



P[bar] 


A/3i[10-3] 


A/32 [10-3] 


A/33[10-3] 


A/34 [10-3] 


A/35 [10-3] 


A/3abm[10-3] 


A/3bw[10-3] 


12 


-0.378 


-2.360 


-2.047 


-3.292 


-1.374 


-7.026 


-4.975 


16 


-0.642 


-3.004 


-2.559 


-4.033 


-1.868 


-8.904 


-6.466 


20 


-0.719 


-3.355 


-2.862 


-4.538 


-2.131 


-10.024 


-7.252 


24 


-0.821 


-3.106 


-2.643 


-4.563 


-2.637 


-10.306 


-7.209 


28 


-0.957 


-3.074 


-2.623 


-4.670 


-3.059 


-10.804 


-7.482 


34.4 


-0.688 


-3.069 


-2.670 


-4.999 


-2.907 


-10.975 


-7.282 



The mean field P-T phase diagram in the case with no anisotropic scattering effect is shown in FiglT] where the 
solid lines are Tab{P) in the impure case and result from the use of the value 1/{2ttt) = 0.13 (mK), while the dashed 
line is that for the bulk liquid (i.e., in clean limit). In obtaining the curves in this figure, we have used the same 
experimental data as in the previous work p3| . Then, a necessary set of Landau parameters was available only at 
three pressure values, 24, 28, and 34.4 (bar) above the bulk polycritical point (PCP). This is why the resulting bulk 
Tab{P) curve has an inessential kink in the figure. The lower (upper) solid line, implying an AB transition curve, 
is obtained without (with) the impurity-induced SC correction given by eq.(18). As one can see in Figl?! the ABM 
region shrinks as the impurity scattering is enhanced. It is important to note that the Tab line obtained with (5/3j 
lies at much higher pressures than the other corresponding curves obtained in clean limit, in the spin-fluctuation 
approach, and in the approach used in the previous subsection. It suggests that the use of this phenonienological 
approach explains better the experimentally determined global phase diagrams of liquid ^He not only in clean limit 
|14| but also in the weakly disordered case, i.e., in aerogel. 

III. APPLICATION TO UNIAXIALLY ANISOTROPIC AEROGEL 

Since the coherence length of superfluid ^He at T = becomes shorter as the pressure increases, can be 
much shorter than the structural correlation length of aerogel ~ 30 — lOOnm at high pressures (see sec. I). In this 
high pressure region, the orientation of scatterers is not random over the scale ~ ^q, and the Cooper pairs can be 
regarded as behaving in scattering potentials with a fixed uniaxial anisotropy. When this local anisotropy induced by 
the aerogel structure is well defined, the role @ of aerogel acting as a "random pinning" of the order parameter field 
(such as the Z-vector in the ABM state) is a minor correction to the free energy, and a mean field analysis assuming 
a fixed anisotropy in the scattering events is expected to well describe true pairing states of the A-like and B-like 
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FIG. 7: The corresponding result to Fig. 6 in the phenomenological approach obtained in terms of l/(27rr) = 0.13 (mK). In 
comparison with the case of the perturbative approach, the shrinkage of the ABM pairing region is much more remarkable. 

phases in aerogel which is locally anisotropic but globally isotropic. [§1 Or, it is possible that an aerogel sample used 
in experiments has accidentally or artificially a global anisotropy. Commonly to these two situations, the mean field 
analysis assuming a small but nonvanishing anisotropy in the impurity-scattering amplitude becomes an appropriate 
description Based on these facts, we apply the expressions of the impurity induced SC corrections to the case of 
an uniaxially anisotropic aerogel by introducing 8] an anisotropy in the scattering amplitude. In performing this, the 
effect of the anisotropy on the SC correction is neglected, because both of them are small effects, and we have only 
to incorporate 60j and 6f3j in the (3j parameters in the analysis performed in Ref.[8|. Thus, only the obtained results 
will be discussed in the remainder of this section. 

In the previous paper Q we found that, in the T-P phase diagram, the ABM pairing region becomes wider in both 
uniaxially stretched and compressed aerogels and, especially, that the polar pairing state appears near Tc{P) as a 
3D phase in the uniaxially stretched case. As is easily expected, the stability of the ABM pairing state is weakened 
by including the impurity induced SC correction SPj and, as a result, the ABM pairing region will shrink. In Fig. 8, 
a result of the uniaxially stretched case based on the perturbative approach is shown, where the parameter values 
{2ttt)^^ = 0.13 (mK) and A = 0.821 are used together with the anisotropy parameter [1| ^« = —0.06. Since the 
polar pairing state is induced by the quadratic term of the GL free energy functional, its region is unaffected by the 
inclusion of df3j in contrast to the shrinkage, indicated by an arrow, of the ABM pairing region at high pressures. 
A similar feature can also be seen in the uniaxially compressed case. Further, although we have shown here only 
the result in the perturbative approach, we have verified that essentially the same result as above is obtained in the 
phenomenological approach. 



FIG. 8; Phase diagram in the perturbative approach for the uniaxially stretched case obtained in terms of (Ittt)'^ — 0.13 (mK), 
A — 0.821, and 5u = —0.06. Since the polar pairing state is determined by the quadratic term in GL free energy functional, its 
region is unaffected by changing the SC correction in f3j. In contrast, the ABM pairing region shrinks due to the inclusion of 



Next we investigate the case with a weaker anisotropy in order to obtain a comparable result witli ^He in a realistic 
aerogel. We choose the anisotropy value 5„ = —0.005 and compare the result in the perturbative approach and that in 
the phenomenological approach with each other. Those results are shown in Fig. 9. Commonly in both the approaches, 
the magnitude of the anisotropy |(5„| is so small that the polar pairing region is too narrow to become visible in the 
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FIG. 9: Phase diagrams in the uniaxially stretched case obtained in the perturbative approach (left) for {2-kt) ^ — 0.15 (mK), 
A — 0.821, and 5u ~ —0.005 and in the phenomenological approach (right) for (27rr)~^ = 0.15 (mK) and 5u = —0.005. Each 
arrow indicates a shift of rAB(-P) due to the inclusion of SPj. Due to its inclusion, the slope of Tab(-P) can easily become 
positive at least in the phenomenological approach. Although the slope in the left figure is not positive, a positive slope can be 
achieved by increasing value for the fixed 5^ even in the perturbative approach. 

T-P phase diagram. The ABM pairing region survives even at lower pressures than the mean field "triangle" (see 
Figs. 3 and 7), and the PCP doesn't exist. At high pressures, the slope dTx^/dP becomes not negative but positive. 
It suggests that the previous result, such as Fig.l in Ref.01, at higher pressures will be improved if (5/3,- is included 
there. An important difference between the two approaches is also seen in Tab{P) at the highest pressure, P = 34.4 
bar, in the figures. Reflecting the contribution of Sf3j, Tab at 34.4 (bar) in the right figure of Fig. 9 following from 
the phenomenological approach is higher than that of the bulk liquid. It implies that the features of the T-P phase 
diagram obtained in the phenomenological approach [l3 | are in better agreement with experimental ones than those 
in other approaches. 

Furthermore, the shaded region of the right diagram in Fig. 9 resembles the experimentally obtained region of the 
A-like phase in Ref.4 where dTAB/dP is positive even at much higher pressures than PCP estimated there |^]. It 
might imply that a global anisotropy in aerogel samples, created through their production processes, needs to be 
taken into account in understanding the experimental phase diagram. That is, it is possible that the window of the 
A-like phase at lower pressures is too narrow to be observed experimentally. Alternatively, at lower pressures where 
Co 3> may be satisfied, a Tc shift due to the quenched disorder [0| outweighs that due to the anisotropy, and the 
scenario in Ref.9 invoking the A-like phase induced by the quenched disorder may be more appropriate. 

IV. CONCLUSION 

We have investigated the impurity- induced strong coupling correction up to the lowest order in l/rT by taking 
account of additional contributions overlooked in the relaxation time approximation and have calculated such contri- 
butions to the quartic term of GL free energy functional using the perturbative approach in a short range repulsive 
interaction and the phenomenological approach. In both approaches, the impurity scattering favors the BW pairing 
state, and the ABM region shrinks just as suggested experimentally The resulting extent of the shrinkage in 

the phenomenological approach [l^] is especially substantial and seems to be consistent with the expectation Q. 

Here, several approximations used in our calculations will be discussed. First, our derivation of the AB transition 
line based on the GL free energy with no sixth order term is not conventional. Although this treatment might have 
a subtle aspect in obtaining the bulk AB transition curve far apart from Tc{P), it is safely valid for the present 
purpose of obtaining the AB transition curve of the liquid '^He in aerogel occurring near Tc. In our treatment on the 
quasiparticle scatterings, angular dependences of the scattering amplitude were neglected. For the scattering events 
in the locally anisotropic aerogels, this approximation of the s-wave scattering may be too simple, and we cannot 
necessarily exclude a possibility that it might have led to different angular dependences significantly affecting the 
impurity-induced SC correction 6Pj. Further, only the lowest order terms in the impurity scattering were kept in our 
calculation by assuming Tr > 1. Although this treatment is safely valid for the present purpose focusing primarily on 
the high pressure region in which the SC effect is more important, in low enough pressures along the Tc(P)-curve or 
close to the quantum critical point, the above inequality is not satisfied, and rather, calculations need to be performed 
in the dirty limit. 

After performing microscopic calculations in the case with an isotropic scattering amplitude according to the two 
approaches, the present theory has been applied to superfluid ^He in aerogel by assuming the aerogel sample to have 
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an uniaxial anisotropy in scattering events over length scales longer than the coherence length ^o- We find that the 
unexpected positive slope of the AB transition curve [1] can be obtained at high pressures, and that the ABM pairing 
region survives even at low pressures with a narrow width. The assumption of an uniaxial anisotropy over large 
scales is realistic from two different points of view. First, as mentioned in sec. I, the relation *C may be satisfied 
in a globally isotropic aerogel under higher pressures, and then, the condensation energy is well approximated Q 
by that in a globally anisotropic case. Second, aerogel samples can naturally or artificially (20l . [2lj | have a global 
anisotropy as a result of their production process. For either of such two situations, the phenomenological approach 
gives more consistent results with experimental facts compared with other models. Consequently, we feel that 
an identification j^, |22| between a superfluid glass phase with the ABM pairing and the A-like phase has obtained an 
additional support. 
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